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The paper deals with the nonequilibrium two-lead Anderson model, considered as an adequate 
description for transport through a d-c biased quantum dot. Using a self-consistent equation-of- 
motion method generalized out of equilibrium, we calculate a fourth-order decoherence rate 7" 
induced by a bias voltage V . This decoherence rate provides a cut-off to the infrared divergences 
of the self-energy showing up in the Kondo regime. At low temperature, the Kondo peak in the 
density of states is split into two peaks pinned at the chemical potential of the two leads. The height 
of these peaks is controlled by . The voltage dependence of the differential conductance exhibits 
a zero-bias peak followed by a broad Coulomb peak at large V, reflecting charge fluctuations inside 
the dot. The low-bias differential conductance is found to be a universal function of the normalized 
bias voltage V/Tk, where Tk is the Kondo temperature. The universal scaling with a single energy 
scale Tk at low bias voltages is also observed for the renormalized decoherence rate j^'/Tk- We 
discuss the effect of 7" on the crossover from strong to weak coupling regime when either the 
temperature or the bias voltage is increased. 
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I. INTRODUCTION 

Over the last ten years, an intense experimental and 
theoretical activity has been developed to study quan- 
tum dots. Due to the presence of strong electronic corre- 
lations in the dot, these mesoscopic systems give rise to 
rich collective phenomena such as the Coulomb blockade 
and the Kondo effect. Their manifestations in transport 
can be studied in a detailed and controlled way in such 
devices as semiconductor-based quantum dots embedded 
in a two-dimensional electron gasi or carbon nanotubes^. 
One of the great interests of these systems is to offer the 
possibility of studying them under nonequilibrium con- 
ditions when either a bias voltage is applied to the leads 
or an electromagnetic field irradiates the device. 

A simple model describing quantum dots is the Ander- 
son model^, in which the dot is represented by a local- 
ized level connected to Fermi seas of conduction electrons 
through tunneling barriers. When the dot is singly oc- 
cupied, it has been shown that the linear conductance 
increases as one lowers the temperature, and eventually 
reaches the unitary limit 2e 2 /h in the case of symmetric 
coupling to the leads. This was predicted in the context 
of quantum dots twenty years ago^ 5 - and was observed 
experimentally about ten years later— ^. 

While in equilibrium most of the properties of the 
Kondo effect are now well understood- 6 thanks to the de- 
velopment of a panel of powerful techniques (e.g. renor- 
malization group, Bethe ansatz, Fermi-liquid theory, 
conformal field theory, density matrix renormalization 
group, slave boson and equation-of-motion approaches), 
most of these techniques fail out of equilibrium. Hence 
there is a huge interest to develop new techniques to 
tackle the problem of the Kondo effect out of equilibrium 



and more generally nonequilibrium effects in strongly- 
correlated electron systems. 

Theoretically, the Kondo effect out of equilibrium has 
been investigated by a variety of techniques developed 
most of the time within the Keldysh formalism: per- 
turbation theory and perturbative renormalization group 
approach^—, slave-boson formulation solved by us- 
ing either mean-fielcU^ or non-crossing approximation^, 
equation-of-motion approache d 15 ' 16 . Exact solutions at 
the Toulouse limit have been proposed^. Other ones 
have extended the Bethe ansatz out of equilibriu m 18 : 19 , 
and in some cases have used the results to construct a 
Landauer-type picture of transport through the quantum 
dot. There have also been important efforts to develop 
numerical techniques such as time-dependent Numerical 
Renormalisation Group 20 ' 21 and imaginary-time theory 
solved by using Quantum Monte Carlo 2 ^. All those ap- 
proaches have only a limited validity of their parameter 
regimes since they mostly describe the properties of the 
system in its ground state, and not in its excited many- 
body states reached when the bias voltage drives a cur- 
rent through the dot. 

In this paper, we develop an equation-of-motion 
(EOM) approach to tackle the nonequilibrium Kondo ef- 
fect in quantum dots and study the decoherence effects 
induced by a bias voltage. The EOM method, though 
conceptually simple, requires some care. The recursive 
application of the Hciscnbcrg equation of motion 2 ^ gen- 
erates an infinite hierarchy of equations, which relate the 
different Green functions of the system. This hierarchy 
has to be truncated by a suitable approximation scheme 
in order to form a closed set of equations. The choice of 
the truncation scheme is crucial in order to treat carefully 
the correlation effects both from the Coulomb interaction 
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and from the dot-lead tunneling. 

The EOM technique was applied to the original An- 
derson model at equilibrium a long time ago^ir— in the 
context of the dilute magnetic alloys. When applying 
the standard approximation based on a truncation of the 
equations of motion at second order in the hybridization 
term t a , it yields results which agree with perturbation 
theory calculations for temperatures above the Kondo 
temperature, Tk- This truncation scheme is usually re- 
ferred to as the Lacroix approximation^. Even though 
the scheme has serious drawbacks at this level of approx- 
imation (underestimation of the Kondo temperature Tk, 
absence of Kondo effect just at the particle-hole symmet- 
ric point), it is acknowledged to provide a valuable basis 
for the description of the Kondo effect both at high and 
low temperatures. The applicability of the Lacroix ap- 
proximation is nicely reported in a recent paper by V. 
Kashchcyevs el a$&. 

In the early nineties, Meir, Wingreen, and Lee2£ under- 
took to apply the EOM method to the study of quantum 
dots out of equilibrium and/or in the presence of a mag- 
netic field. They used a simplified version of the Lacroix 
approximation, which fails to account for the finite de- 
coherence rates induced by bias voltage and/or magnetic 
field. Meir et at proposed to introduce them hcuris- 
tically by making use of the Fermi golden rule. They 
obtained interesting results for the bias voltage depen- 
dence of the differential conductance showing a zero-bias 
anomaly, but the weakness of the approach is that it 
does not constitute an unified and consistent frame for 
the treatment of the Kondo effect in the presence of the 
dccoherence effects induced out of equilibrium. 

There have been recent attempts to use an approxima- 
tion which truncates the equations of motion at higher or- 
der in t u 15 i 16 i 30 . Their authors claimed to improve quan- 
titatively at equilibrium the Kondo temperature and the 
density of states around the Fermi energy and have been 
able to investigate some nonequilibrium issues. How- 
ever, there is need to clarify the decoherence effects in 
the framework of the EOM method. 

The organization of the paper is the following: 

- In Sec|nl we outline the EOM formalism and the 
main steps of the proposed approximation based on a 
truncation of the equations of motion at the fourth or- 
der in t„. An analytical expression of the retarded 
Green function in the dot is derived involving expecta- 
tion values which are determined self-consistently. The 
details of the calculations are presented in |Appendix I| 
and | Appendix II} 

- An analytical study of this Green function is pre- 
sented in Sec lIIIl Namely, we deduce the renormalization 
effects and the transition rates involved at the second or- 
der in t a for the different regimes of the Anderson model. 
Special care is given to the singly-occupied dot regime, for 
which the vanishing of one of the transition rates leads to 
the low-energy logarithmic divergence of the self-energy 
of the dot Green function, yielding a Kondo resonance 
peak in the density of states. We show how the approx- 



imation scheme allows one to derive a decoherence rate 
out of equilibrium which becomes finite as soon as a bias 
voltage is applied. This decoherence rate provides a cut- 
off to the logarithmic divergence present at equilibrium 
in the Kondo regime. We show how our approximation 
scheme improves the result for the Kondo temperature 
Tk upon earlier predictions. 

- We present in Sec lIVI our numerical results - both at 
equilibrium and out of equilibrium - for the density of 
states, the (linear and differential) conductance and the 
bias-induced decoherence rate. A self-consistent treat- 
ment is required in order to determine the expectation 
values involved in the Green function of the dot. At 
equilibrium, the density of states in the particle-hole sym- 
metric case shows a three-peak structure at low temper- 
ature with a Kondo resonance peak in the local moment 
regime. This result constitutes an advantage of our ap- 
proximation scheme compared to the Lacroix approxi- 
mation. The unitary limit for the linear conductance 
G = 2e 2 /h is analytically recovered at zero temperature 
in the particle-hole symmetric case when the dot is sym- 
metrically coupled to the leads. Numerically, one notices 
a slight underestimation of G due to the numerical accu- 
racy of the self-consistent treatment. 

Out of equilibrium, the spectral function shows a split- 
ting of the Kondo peak into two peaks pinned at the 
chemical potentials of the two leads. The height of the 
peaks diminishes when the bias voltage increases, mean- 
ing that the Kondo effect is destroyed by decoherence 
induced out of equilibrium. This influences the evolution 
of the differential conductance as a function of bias volt- 
age, which shows a zero-bias anomaly followed by a broad 
Coulomb peak. At low bias voltage, we check that the 
differential conductance follows a universal scaling law 
which depends on a single energy scale, Tk- Comparison 
is made with results obtained by recent numerical tech- 
niques for nonequilibrium such as time-dependent NRG^I 
and imaginary time theory solved by QM C 22 ' 31 . 

Finally, we compute the bias voltage dependence of the 
decoherence rate and discuss the crossover from strong 
coupling to weak coupling regime depending on the com- 
parison between the decoherence rate and a characteris- 
tic energy scale T* . We show the existence of a crossover 
between the strong coupling and weak coupling regime 
when either the temperature or the bias voltage is raised, 
as predicted by other methods. 



II. EQUATION-OF-MOTION FORMALISM 

We model the quantum dot connected to the two leads 
by the single-level (spin-1/2) Anderson impurity hamil- 
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tonian 



H 

Hlead 
Hdot 
Hi n t, 



Hlead + Hdot + Hint, 



(1) 



Sang + Un^ni, 

a 



where c aka (c ak a) is the creation (annihilation) operator 
of an electron of momentum k and spin a(— ±1) in the 
a(= L, R) lead (with energy e ak — E k — Ma); Ma is the 
chemical potential in the a lead; ft(fa) is the creation 
(annihilation) operator of an electron of spin a in the 
quantum dot (with energy e a — s c i + uB/2 when the 
Zeeman splitting B is included); n a = ft fa is the num- 
ber operator for electrons of spin a in the dot; U is the 
Coulomb interaction between two electrons of opposite 
spin in the dot; and t atT is the tunneling matrix element 
between the state \ka) in the a lead, and the state \a) 
in the dot. For simplicity, we assume t aa to be real and 
fc-independent. 

When a bias voltage is applied to the leads (eV = 
Ml — Mi?) j the system is driven out of equilibrium and a 
current is induced through the quantum dot. The current 
/ for the Anderson model is expressed by the general- 
ized Landauer formula^ accounting for the interactions 
among electrons 



, r/ J(T (e)rfl (T (e) . L R 
r — TTTT — rT w e )-/p(e) P»- £ ), 

W 1 La{£) + I Ra{£) 



(2) 



where 



• W is the half-bandwidth of the conduction electron 
band in the leads, 

• T a(T (e) is the tunneling rate of the spin a dot 
electron at energy e into the lead a, defined as 
Y aa {e) = nJ2 k t 2 a J(e - e ak ) = nt 2 ary p° a (e) with 
Pa(e) the unrenormalized density of states at en- 
ergy e in the lead a, 

• fp(e) — {cxp \j3(e — p a )] + is the Fermi-Dirac 
distribution function in the a lead, 

• Po-^s)-, the local density of states for spin a in the 
dot, can be expressed in terms of the retarded elec- 
tron Green function in the dot Q T a {e) according to 
p a (e) = -l/nlmg r a (e). 

As pointed out in RefHU, Eq.© is valid provided that 
the tunneling couplings for both leads Ti ja (e) and r^ cr (e) 
differ only by a constant multiplicative factor. 

The task is to compute the retarded electron 
Green function in the dot defined as <?£(w) = 
-i / °° dte^({f a (t), /t (0)}) where C - w + iS (6 0+). 
In order to simplify the notations in the rest of the paper, 
the imaginary part id going alongside u> will be implicit, 
while the summation over k implies summation over both 
a and k. Hence we write in a shorthand notation 



Yj ta °^12 ] 



Using the Zubarev notation^ for the retarded Green 
functions involving fermionic operators A and B 



((AB)) 



-i lim 

5^0+ .,o 



dte i{u+i6)t ({A(t),B(0)}), 



. (3) 

which can be integrated by parts, and using the Heisen- 
berg equation of motion, one can show the following re- 
lation 



uj((A,B)) = ({A,B}) + (([A,H],B)). 



(4) 



This allows us to derive a flow of equations for the dot 
Green function 7 a (uj) = ((fa, ft))- We also adopt a sim- 
pler notation in the following derivations by changing 

((A, ft)) ^ ((A)). 



Applying Eq. we find the first equations of motion 
(w -e*)((U)) = 1 + J2^(( c k,))+U((n s f a )),(5) 



(uj - e k )((c ka )) = t a ((f a )). 



(6) 



Combining Eqs. (J5J [6]) yields 

[w - e a - E°(w)]((/ ff )} = 1 + U((n s f a )), (7) 



where S°(w) = 



Eqs. (|5|6p are referred to as 



the first-generation equations of motion in the hierarchy. 
They govern the evolution of the Green functions formed 
by a single operator (i.e. ((fa)) and ((c krJ ))). 

Throughout this paper, we assume that the half- 
bandwidth W is much larger than all the other energy 
scales, so that the band edge effect does not affect the 
local density of states in the dot p a {u)- In this case, the 
properties of the system at low temperatures do not de- 
pend on the exact value of W since only states around 
the Fermi level contribute, justifying the consideration 
of the wide-band limits W — > oo. Within this limit, 
the non-interacting self-energy can be approximated by 
— — iT a , where T a (= + Tj^ a ) is a constant, 
independent of energy 

Interesting dynamics comes from ((n a fa))\ its equation 
of motion is given by 

[w — e a - U]((n s f a )) = (n s ) + ^ \ta((n s c k a)) 

k 

+ta((.ftckafa)) -ta((clfafa))} - (&) 

The Green functions appearing on the right-hand side 
of Eq. © have their evolution governed by the following 
equations 



a=L,R k 
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(u>k:aa - U) ((4_/^/ CT )) 



tAMa)) + ^[<</ < W W >) - ((c{, s f & C k<T )) 
fc' 

(flc kS ) +t s ((n s f a )) + Y,[t«(UUk S c k >«)) - ts((4<teW) 
fc' 



(9a) 
(9b) 
(9c) 



r 



where we write in a shorthand notation 



£b 



with {a/3 • • • , ab ■ ■ ■ } being any set of parameters within 
k's and cr's. We have for instance: uj :k = u—£ k , ^k- = u+ 

£k, ^a-.ka = U + £a-£k~£a, &nduJ k . aS = U + E k - E a - Eg . 

Eqs. © generate three new Green functions on their 
right-hand side. Generally one can divide the expressions 
of the latter Green functions into two parts 



((fUkaCk'a)) 
({cbfvCk'v}) 



{{c\ s Ck>gfa)) 



(ftc k t)((c k , a )) + ((/IctoCfc^UlOa) 

(<LU)(M) + «4/^>Uiob) 

(cLc k , s )((f„)) + ((cl-c k ^f„)) c {10c) 



where the first part is obtained by decoupling pairs of 
same-spin operators and the second part ((• • -)) c de- 
fines connected Green functions in the spirit of cumu- 
lant expansion. It is often assumed that these connected 
Green functions arc negligible. This assumption has 
been broadly applied, usually referred to as the Lacroix 
approximation^ . It turns out that within this approx- 
imation, the calculation of ((fa)) is exact at the second 
order in t a , while it picks up some of the fourth-order 
contributions. At zero temperature, this approximation 
leads to logarithmic singularities in the density of states 
of the dot at the chemical potential, even in the pres- 
ence of an external magnetic field B or a bias voltage V. 
These divergences arc unphysical since one expects the 
logarithmic singularities to be washed out by decoherence 
effects introduced by either B or V. 



In this work, we propose to go beyond the Lacroix ap- 
proximation, and consider higher hierarchy in the equa- 
tions of motion. The interest of the approach is to 
account for the decoherence effects introduced by ei- 
ther nonequilibrium or the presence of a magnetic field. 
The detailed derivation of the higher-hierarchy equa- 
tions of motion and the decoupling scheme are given in 
Appendix I In the approximation scheme we propose, 



after having expanded the equations of motion to order 
t a , we decouple pairs of same-spin lead electron operators 



,t 

'fc'er 



Cfccr)) and pairs of same-spin dot-lead electron 



(e.; 

operators (e.g. (c\ a .f a ) and (/]c fecr )). This decoupling 
has to be done carefully in order to avoid double count- 
ing. The equations of motion of the three functions on 
the left- hand- side of Eqs. (JSJ are given by Eqs. (|A.12[) . 
which are exact up to order t a . 

Combining Eqs. ([71151 1 A. 12j) yields a rather complex ex- 
pression for ((f a )) since Eqs. (|A.12[) couple to each other 
in an integral way. In this paper we will limit ourselves to 
a simple expression for the Green function ((fa)) that ne- 
glects contributions generated through integral coupling 
of the equations of motion. This is motivated by the fact 
that the contributions we keep effectively lead to a second 
order term after resummation, as shown later on. The 
integral terms we neglect, on the contrary, are at least 
of fourth order in t a and we believe they are not rele- 
vant for our study of the decoherence effects at nonzero 
bias and/or temperature. Thus, we are able to reduce 
Eqs. (|ATT2|) to 



[U3..k - £i CT (u; : fc)] ((n S C k a)) 
[u S :ka - Z 2 a{oJ:k)] ((flckafa)) 

[U k ..aa - U -%aM]((c{jafa)) 

[u. k - Yli a (uj :k )] ((n a c k a)) 



ta((n S fa)) + S 5 a(w : fc)((naCfc CT )), (11a) 

(fUka) + t s ((n s fa)) - £ [fe/Sfc + (fUka^D S ;kk> ]T ((fa)^) 
fc' fc" 

(Aafa) - t S ((n S fa)) + £ [WfcV + (4^)4^:* £ Fk» k] «/<r)>(Hc) 
fc' fc" 

-(PaC k a) + £ [taf° k , k + (/Jcfa^iW £ «/„)) 



+ Z 5 a{U:k)((naC k a))- 



(lid) 



where we define f k , k = (c\ s c k >a), 
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S lff (a; :fe ) = £4(^+^4*). 

k' 

fc'fc" 

+ £(4^:fctf+4^U)> (13) 
fc' 

k'k" 

+£(4^4- +4^), (14) 



E 5o -(w :fc ) = a; :fc t g ^fc' :a -fc(4' g /^) + D cf.kk' (/]cfcv) 
fc' 

+ £ 4(D fc , CTfc / fc V - Da-kk' fk"k') (15) 



(12) In the last equation Eq. (fT6)). the sign in front of U is the 

same as the sign in front of e a in ui a ^—:ab--- ■ Thus we 

have for instance: D a .. kk , = -Uu~^ kkl (uJa-.kk' + Uy 1 = 

-U (w + e a - e k - ek')' 1 (u> + e a - e k - e k > + U)~ l . 

Notice that we keep heuristically a fourth-order term 

on the right-hand side of Eq. (jlld|) — explicitly the 

term (fUka)Y^k'k/'( t l D ^'fk''k')((f^)) ~ in order 
to respect the unitarity condition typical of a Fermi 

liquid Im[Cf£] -1 = T a . This is explained in detail in 

Sec lIV Bl We emphasize that in principle this term 

shall be recovered by properly including fourth-order 

contributions generated through integral coupling of the 

equations of motion (|A.12j) . 

Combining Eqs. ([JJ [5J [11]) yields the following expres- 
sion for the Green function in the dot 



fc'fc'' 



and 



D a p... :ab ... = -Uu> a l.. ab ...(uj a { 3 ..... ab ...±U) 1 . (16) 



U2„{u) - {n s ) + U a (u>) 



(17) 



where we define the functions 
u la {uj) = UJ.a- - E°(w), 



1 

u 



ti 



UJ- k - E 6(7 (w : fc) cj-. k(T - S 2CT (w : fc) -UJk-.ag- + U + E 3<T (w fe: ) 



^6cr(^:k) = T, lcr (uj. k ) + 



E 5g (u :fc )E 5cr (u :fc ) 
L0 :k - Ei S (w :fc ) 

ta(fUk<r) 



-£ 



* cr ( C fca/ <7 



fc ^ 



5:fc CT - S 2CT (o; :fe ) Y -^^ff + C/ + E 3a (w fe: ) Y ^ ~ E 6 CT (W:fc)][W:fe - Elff(w.fc)] 



E 



t CT E 5(f (w :fe )(/J;cfc (T ) 



(18) 
(19) 

(20) 
(21) 



s,(«) = -£ 
fcfc' 

+£ 



— I +E i — 



fcfc' - S 2<x(^:fc) fcfc/ 

^ff /fc'fc + tcr(f^C kcr )tlD a:kk > Y,k" fk"k'] S 5<r(^:fc 



-Wfc :CTff + £/ + E 3o .(w fc: ) 



fcfc' 



[u) :k - E 6(T (u; :fc )][w :fe - Ei ff (w : fe)] 

I 



(22) 



In order to close the problem, one needs to append to 
Eq. flTTl) the closure equations, which enables one to deter- 
mine the expectation values showing up in the expression 
of Q r a (u>). The calculations of the expectation values are 
presented in |Appcndix II[ while the full self-consistent 
treatment is explained in Sec lIV A"l G^{u>) given by Eq. 
(|17[) respects charge conjugation symmetry, as proved in 
|Appendix III| 



regimes of the Anderson model, as can be derived from 
the results obtained in the previous section. Special care 
is given to the singly-occupied dot regime where many- 
body effects can give rise to Kondo physics. We ana- 
lyze in detail the nonequilibrium situation in the latter 
regime, and show how the EOM method provides a pow- 
erful frame to describe the decoherence effects induced 
when a bias voltage is applied to the leads. 



A. Renormalization effects and transition rates 

III. ANALYTICAL RESULTS 

In the presence of the Coulomb interaction U and the 
In this section, we discuss some aspects of the behavior dot- lead tunneling coupling T acr , the bare parameters of 
of the system in and out of equilibrium for the different the Anderson model get rcnormalizcd according to (for 
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Tier 




f 21 




Empty dot 
(e CT - fie > r) 


r 


r 


r 





Kondo regime 
(e CT + U - fi a , fJ-a - £o > r) 


r 





2r 


r 


Doubly-occupied dot 

Oa - e CT - u » r) 


r 


r 


r 





Mixed valence regime 
(Min{£ CT -fi a ,/i a -t/} ~T) 


r 


r 


r 






TABLE I: Transition rates 7P 2 ' = — ImE; at the second order in t CT and at zero temperature, for the different regimes of the 

Anderson model obtained by the EOM approach. Notice that, in the Kondo regime, 73 = yields low-energy logarithmic 

(2) 

divergence of the self-energy of the dot Green function, responsible for the Kondo effect. In the latter regime, 7g CT ^ 0, which 
brings on an additional divergence arising from Eq. (|25|) . 



zero temperature) 

4 - E if- U^feTSl _„.„ ), (23.) 

The above results are obtained from Eq.(fT7| up to sec- 
ond order in t a and by taking f£ k , = fp(£k)&kk' ■ In the 
mixed valence regime (Min{e cr — (i a , fi a — e a — U} « T), 
the renormalization of the bare level energy is consistent 
with the prediction of the scaling theory^^ as pointed 
out in the previous EOM studies 2 ^ As expected, the 
renormalization effects are small around the particle- 
hole symmetric case (e a = —U/2). In the large U/\e a \ 
limit, the renormalization effects are very important, as 
it is the case for quantum dots coupled to ferromagnetic 
leads^. 

Interestingly, these renormalizations are consistent 
with the shift of the pole of the Green functions 
({fickafa)) and ((c\ a f & f a )). For instance from Eq. 
(|llb[) . the pole of {{fgCkofo)) (with respect to u : k) is 
shifted to, 

£a - £<f + Re S 2(T (£* - £%) = E* - £%. 

The shift of these poles can have important consequences 
on the splitting of the Kondo resonance peak when a 
magnetic field is applied. These corrections are neglected 
in the Lacroix approximation. 

The imaginary part of the corresponding self-energies 
evaluated at the pole of the Green functions (e.g. at 
iu-k = £* a — £% for ((fgCkgfcr))) defines the transition rate 
from the state f a \GS)) to the excited state f a c\ a f s \GS), 
where the ground state is denoted by \GS). Within sec- 
ond order in t a and taking into account the renormal- 
ization of the dot level energies, the transition rates are 



given by 

7ict = -Im£i CT (0) = 2r ffj (24a) 
4 2 J = -Im£ 2CT (4-4) 

= E E r -[W£(<)+/i(<+tnK24b) 

ol—L.R a 

7$ = -ImIl 3a (e* a + 4 + U*) 

= E E r -[ 1 + ^(4)-/F(e; + ^)],(24c) 

a=L,R a 

7 £ 5 = -ImE 5(r (0) 

= 2 y, r aa imo-m< + u*)}, (24d) 

a=L,R 

where T a = r ' i, a + Tfi a using the notation defined in 
Sec[IIl The values of these second-order transition rates 
in the case of spin-independent tunneling (Ef- = = 
r/2) are reported in Table U for the different regimes of 

the Anderson model at zero temperature. One can note 

(2) 

that, in the wide-band limit, the value of 7^ does not 
depend on the occupancy in the dot. In contrast, the 
other transition rates take different values depending on 
the regimes considered. One can distinguish four regimes: 

a) In both the empty and doubly-occupied dot 

(2) n j (2) (2) a (2) „ 

regimes, 7^ ~ and 7^, 7^ ^ E. As 7^ = 
in these two regimes, the third term of ^ a {ui) and 

U a (u}) vanishes (cf. Eqs. I|21l22p ). On the other 

(2) (2) 

side, the finite values of 72(T and 7^ provide a 
cut-off to the integrals involved in the calculation 
of the remaining terms, thereby preventing them 
from diverging at low energy. As a result, the elec- 
tron density of states in the dot does not show any 
resonance peak but only two broad peaks located at 
the positions of the renormalizcd dot level energies. 
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b) In the mixed valence regime (take for instance e a — 
fJ-a ~ T), the renormalization effects push the dot 
level energies above the chemical potential, hence 
the transition rates are identical to those found in 
the two regimes of a) . Our numerical results for the 
density of states are in better agreement with the 
exact numerical renormalization group result than 
those found in the Lacroix approximation or the 
non-crossing approximation, for which a spurious 
peak may ap pea r at the Fermi level, as it has been 
shown in Reffla. 

c) The singly-occupied dot (Kondo) regime is the 
most interesting since one of the transition rates 

(2) 

72 CT vanishes. This gives rise to a logarithmical 
divergence at low energy of the integral involved 
in the calculation of the first term of S<t(w) and 
H a (uj) in Eqs. (|21I22[) . Another divergence comes 
from the calculation of the third term of S CT (u;) and 

(2) 

H a (uj), which no longer vanishes as is now fi- 
nite. The integrand of those terms has a structure 
like 

S 5 ,r(W:fe) 



[Uj-k - 2^ 6a (uj. k )}[uj. k - Si ff (w :fe )] 



r 



i 



i 



UJ-.k U>:k + 2iT 



o 



(25) 



+ *76oa and 

7^ = and 



These terms have two poles at uj- k 

n i ■ ( 2 ) - + i ( 2 ) (2) 
uj :k = + with = 7 y 

(2) (2) (2) 

ibba = ns +75<r = ^r, respectively. The values of 
the imaginary part of these two poles are reported 
in Table U for the Kondo regime. Since 7^ = 0, 
the first term in Eq. (|2"5"|) gives rise to an additional 
logarithmical divergent self-energy term at low en- 
ergy. The presence of these two logarithmical di- 
vergences mentioned above is solely responsible for 
the formation of the Kondo resonance peak in the 
electron density of states in the dot. In Sec lIII B"3l 
we will analytically estimate the Kondo tempera- 
ture from the consequences of these divergent self- 
energy terms. 



B. Case of the Kondo regime 

The Kondo regime is particularly interesting because 
some logarithmical divergent terms (Kondo singularities) 
survive even after introducing second-order self-energy 

corrections £) a :{i = 1,2,3,5), as discussed before. We 
focus in more detail on this regime and show how fourth- 
order corrections in t a smear the Kondo singularities 
when the system is driven out of equilibrium. 

1. Compact expression for the electron Green function in 
the dot 

To facilitate the understanding of Eq. (TIT)) , it is in- 
structive to put the expression of G 7 a (w) in the Kondo 



regime in a more compact way in order to better identify 
the terms bringing about Kondo singularities. After in- 
tegrating over k and using the closure equations for the 
expectation values (see | Appendix Ilj ), we can express the 
functions E^w) and ILj.(cj) appearing in <?J(w) as 

S a (oj) = -7 + oQ ff ( w )-ff ff roP ff ( w ), (26) 
U a (u;) = roP,(«). (27) 

where the functional J 7 acting on any function A" ct (oj) 
is defined as 



±XA -oj :(T ,f* +U* -2iT 



r 



-XAu + i^) +XAu + 2iT 



(28) 



where 'cr*' in u>g* :a * and La- a - s - indicates that the dot 
level energies are renormalizcd. 

The derivation of P a {ui) and Q a (u>) is given in 
|Appcndix II| while the full self-consistent treatment is 
discussed in Scc llV Al We report here the result obtained 
for P(t(u)) and Q a (cu) 



T O */ 



a=L,R 



uj — e + id 



(29) 



a=L.R 



uj — e + id 



where Q^{s) is the advanced dot Green function. 

One can see from Eqs. (|17I26I27I28[) that the expression 
of G„(w) contains four terms Q s (w ff * :(T * ), P ff (w ff » :(T *), 
Qcr((jj) and P a {ui) which give rise to low-energy Kondo 
singularities when only second-order transition rates are 
considered. We will see in the next section how the 
nonequilibrium situation cures these divergences by in- 
troducing finite transition rates coming from fourth-order 
contributions in t a which provide a cut-off energy to the 
divergent integral terms. 



2. Decoherence rates induced out of equilibrium, 

We calculate explicitly the fourth-order transition 
rates (decoherence rates) by expanding the equations of 
motion to sixth order in t a , followed by the usual trunca- 
tion. The derivation is long but straightforward and we 
present only the results for the fourth-order decoherence 
rates in the Kondo regime, namely, 7^ and 75*^ 
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-Y (4) 
72 a 



(4) 
/oaer 



E E 

a,P=L,R a,a' 

E E 

OL.fi=L,R a,a' 



^- ao^- (3c 



ds(l - /£(e))/£(e -e a + e al )V[D a {e)% 



de{\ - ff(e))f F (e - s a + e a ,)V[D a (e)% 



(31a) 
(31b) 



where 



1 



1 



e — e a + iS e — e a — U + i5 
In the limit V = \pl — Pr\ "C Minde^ — fi eq \ , £d — fJ, e q + U} (with p, eq = (pit, — /xr) /2), and at zero temperature 

T aa' T fi/jft 



(32) 



(4) 
72a 



(4) 



-(/U/3 ~ Ma + £<J> - £ a ")Q{Hi3 - p a + £ a > - e a n)D a >{yL a )D a n{np) 



E E 

a,/3=L,Ro-' -a-" 

^ ^2 E (VP - Pa + £<r' - e a »)®{np - p a + £ a > ~ £a")P° a P°J a a' ,/3a" Jpa'' 
a,P=L,Ro\a" 

E E 

a,p=L,R <y' 

— E _ ^ a £a ' ~ £ cr")©(M/9 — Pa + £a' — £a")PaP°pJaa\Pa"Jl3a",ao 

a,/3=L,R <r' ,<," 



(33a) 



■(flp - p a + £ a > - e a >i)Q(p.p - p a + £„' - £<r» )D <7 > (pa)D a » (flp) 



(33b) 



r 



where 0(x) is the Heaviside step function and V denotes 
the principal value of a function. Eqs. (|33p are expressed 
in terms of the Kondo exchange coupling^ Jaajia' = 
2t atT tpa-i D a (jjLa). In the absence of magnetic field, both 



(4) 



(4) 
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(4) 



decohcrence rates are equal 72 

The expressions of these two decoherence rates are the 
main result of this section. Although for , summation 
is only over opposite spins, both of them involve at least 
one spin-flip process. At zero temperature, these deco- 
herence rates are finite as soon as a bias voltage and/or 
a Zecman splitting is introduced. The finite values of 
these decoherence rates provide a cut-off to the diver- 
gent integral terms of the Green function and smear the 
Kondo singularities. Note that 7t5 is slightly different 
from the heuristical result of Refl29l obtained from the 
Fermi golden rule, because here both spins contribute to 
the rate. Our overall result for the decoherence effect 
is consistent with those found using a real-time diagram- 
matic technique^ and the non-crossing approximation^, 
although in the latter case the decoherence rate was not 
calculated explicitly. 



3. Kondo temperature 

At equilibrium and at zero temperature, the Kondo 
scale Tk (ks = 1) can be roughly estimated from the 
zero of the real part of the denominator of G^(u>) in 
Eq. (jTT)) located near the chemical potential^. Consider- 



ing the case of zero magnetic field and spin- independent 
couplings To- — r ff = r/2 in the wide-band limit, Th- 
reads 



T K * [2r(2 £0 + Uf + 8I-] ^ exp { 47r£ f r ; +L/) } 



(34) 



where £0 = £d~Peq- Tk is independent of W, as expected 
since the high-energy scale is now regulated by U. 

We now compare our result for Tk in Eq.(f34"|) with 
that obtained within the Lacroix approximation [2eo + 
U]cxp [2-Ke (e + U)/TU]. 

First, the Lacroix result for Tk is improved by an ex- 
ponential factor 4/3, in better agreement with Haldanc's 
prediction^ (UT/4) 1/2 exp[n£ (£ + U)/TU}. This is 
due to the presence of an additional logarithmical di- 
vergent term in the self-energy given by S CT (£). To our 
best knowledge, this contribution, coming from a fourth- 
order self-energy, was first found by Dworin^l and was 
attributed to a finite lifetime mechanism of the localized 
electron. It was lately reproduced^ in the infinite U 
limit. 

Secondly, at the particle-hole symmetric point (2eo + 
U = 0), the proposed approximation cures the aforemen- 
tioned pathology of the Lacroix approximation for which 
Tk vanishes, as will be further explained in Sec II V Bl 
The expression of Tk at that point is given by 



T K ^ 2r exp 



27T£ 



3r 



(35) 
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IV. NUMERICAL RESULTS 

We present our numerical results in and out of equilib- 
rium, and discuss the evolution of the density of states 
as well as transport quantities. We consider a quantum 
dot connected symmetrically to the two leads with spin- 
independent tunneling couplings (Ylo = Yr s = Yls = 
T m = T/4), and take a large ratio U/Y < W/Y = 20 in 
order to be in the wide-band limit. For illustrative pur- 
poses, we choose to present the results at the particle- 
hole symmetric point (e^ = —{7/2), which turns out to 
be particularly well described by our method, in con- 
trast with the other EOM approaches developed so far. 
Finally, we limit the study to the case of zero magnetic 
field in order to concentrate on the nonequilibrium effects 
brought by the application of a bias voltage. The Green 
function G"^{oj) given by Eq. (jTT)) is solved in a fully self- 
consistent way (cf. Sec II V AT) . In the Kondo regime, an 
important energy scale is provided by the Kondo temper- 
ature which needs to be properly defined. We will not 
use the approximate expression for Tk given by Eq. (|34[) 
but rather calibrate it numerically from the temperature 
dependence of the zero-bias conductance 



dl_ 

dV 



T=T K ;V=0 



-2^0, 



(36) 



where Gq(= 2e 2 /h) is the zero-bias conductance at zero 
temperature. 



As a result, in equilibrium, Eq. (fl7|) ends up being an in- 
tegral equation with respect to G r a {u) that can be solved 
self-consistcntly. 

However, out of equilibrium, the above relationship 
between the different Green functions no longer holds, 
and one cannot compute the expectation values from the 
spectral theorem. An alternative is to work within the 
Kcldysh formalism. The details of the calculations of 
the expectation values (f^Cka) and (ct a Ch'a) within our 
EOM approach, and of the related integrals P a (u)) and 
Q a (y>) through which these expectation values contribute 
to Eq. |[T7]). are presented in |Appcndix II| (cf. Eqs. IB.41 
IB.5[) . In the wide-band limit, it turns out that even out 
of equilibrium, the integrals P a {^) and Q a (uj) keep the 
same structure as in equilibrium, and depend only on the 
retarded Green function without requiring any knowledge 
of the lesser Green function. 

As far as the occupation number in the dot (n a ) = 
(/ct/ct) i s concerned, the calculation is rather more com- 
plicated out of equilibrium since the simplification which 
takes place before for the calculation of Pa(u) and Qct(cj) 
does not occur, and one needs to know the lesser Green 
function Q< (cj) in order to derive (n a ) by the use of 



duj 
2^ 



(39) 



To find Q<(lu), we use the Dyson equation written in 
the Keldysh formalism G„{oj) = <7£(w)£<(u;)C/£(u;) and 
express the lesser self-energy 5v (w) via the Ng ansata^ 



A. Self-consistency 

The dot Green function given by Eq. (|17[) shows an 
explicit dependence on the expectation values (/Jcfe CT ), 



J 



<'k> 



(denoted by f k t k previously) and 



{flfcr)- What matters then is to compute these expecta- 
tion values in order to properly define the self-consistency 
scheme. In general (in both equilibrium and nonequilib- 
rium situations), the expectation values (as for instance 
(/Jcfeo-)) can be expressed in terms of the related lesser 
Green function 



(37) 



In equilibrium, the relationship Q ka a (oj) 



-/fMI^L^M - Gfo^iu))] holds, relating the lesser 
to the retarded and advanced Green functions, Q r k(J a (uj) 
and Q£ a (uS) respectively. The expectation value is then 
given by 



(fa c ko 



<k)f F (w)lmQ r ka (w). 



(38) 



This relationship is nothing else but the spectral theo- 
rem which expresses the expectation value in terms of a 
functional of the corresponding retarded Green function. 



-2/ 



E ■ 

*=L,R 



-/£(u,)Im£», 



where 5v(w) = w — e a — [^(u)] -1 is the retarded self- 
energy. This ansatz is based on an extrapolation from 
both the non-interacting limit out of equilibrium and the 
interacting limit in equilibrium. Thanks to this ansatz, 
the calculation of (n a ) can be performed from the knowl- 
edge of G^(i-tj) only. Let us also mention that many re- 
sults can be obtained at the particle-hole symmetric point 
(also out of equilibrium), where the occupation number 
is identically 1/2. 

Therefore, all the expectation values relevant to the 
calculations can be expressed in terms of Q^(uj), and the 
self-consistent scheme is straightforward. Eq. (fTT|) ends 
up being again a complex integral equation with respect 
to Q' a (Lo), exactly as in the equilibrium situation except 
that now the different chemical potentials of the two leads 
have to be entered explicitly. We emphasize that this 
constitutes a huge simplification in the technique that 
renders the approach developed in SeclTl] tractable in a 
self-consistent scheme even out of equilibrium. 



B. In equilibrium 

We compute the density of states in the dot p a {^) = 
-l/irlmQ'^Lo) at equilibrium using our EOM approach. 
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FIG. 1: (Color online) Equilibrium density of states in the 
particle-hole symmetric case at T/r = 1CT 3 for different val- 
ues of the parameter U (the chemical potential of the lead fi eq 
is taken equal to 0) . The density of states for large U shows a 
three-peak structure with two broad side peaks and a narrow 
Kondo resonance peak centered at the Fermi level. 




FIG. 2: Linear conductance as a function of dot level energy 
Ed, for U/T — 4 and at different temperatures. When the 
temperature is lowered, the conductance is enhanced in the 
singly-occupied regime — ea/T £ [0, 4], and eventually reaches 
the maximum conductance 2e 2 /h for a single channel at zero 
temperature. The conductance does not reach this limit here 
because of the numerical accuracy of the self-consistent treat- 
ment. 



two broad peaks are centered at the renormalized energy 
levels; their position, intensity and amplitude agree quan- 
titatively with the NRG result^. The Kondo resonance 
peak is pinned at the Fermi level of the leads. 

The fact that our EOM scheme correctly describes the 
particle-hole symmetric case is one of the successes of 
the method. This can be understood by the fact that in 
the previous EOM approaches, for the Kondo regime, 
there is an exact cancellation of the divergent terms 

Q*(U*:*)-QA-U:* S + U) = Qs(u)-Qa{-U) = 0. This 

feature is cured in our EOM approach since the function 
Qs{— uj :a -a + U) in Eq. (|26[) acquires a finite transition rate 
2T and is therefore smeared out. Therefore, the cancel- 
lation does not occur any longer, and we are left with a 
divergence in the self-energy at the origin of the forma- 
tion of the Kondo resonance peak. Through the same ar- 
gument, our approach is shown in Sec lIIIB~3l to improve 
the prediction made previously by the Lacroix approx- 
imation for the Kondo temperature in the particle-hole 
symmetric case. 

Moreover, the density of states at the Fermi level is 
found to be p a (n e q) = 2/ttT in agreement with the Fermi 
liquid property at zero temperature and hence respecting 
the unitarity condition. This can be explained as follows: 
at zero temperature, the functions Pa(uj) and Qa{u>) di- 
verge logarithmically asw-> (i eq , 



n 



^( Me9 )ln| W - Me9 |+0(l), 



Q a {u) = -^[l + *T CT S> e J]ln|cj-/i eg |+0(l). 

We find that the inverse of the imaginary part of <?J(o;) 
(cf.Eqs. (fT7l26l) 1 is Im^]" 1 ^) = T/2, as expected 
from the Fermi liquid theory. In the particle-hole sym- 
metric case, we find R,cQ' a {p, eq ) ~ 0. Combining these 
two results leads to p a (^, eq ) = 2/7rr as observed in FigUJ 
Inserting the value of p<j(p eq ) into Eq. ([2]) allows one to 
find the current at small bias voltages and from there 
the linear conductance G = dI/dV\v=o- When the dot 
is symmetrically coupled to the two leads, the unitary 
limit G = 2e 2 /h is recovered at zero temperature. The 
numerical results for G as a function of the dot level £d 
are shown in Fig. [3J As can be noticed, the method 
underestimates G and the unitary limit is not exactly re- 
covered at Ed = — U/2 because of the numerical accuracy 
of the self-consistent treatment. 



FigfT] reports the result for the density of states at equi- 
librium and T/r = 10 -3 for different values of the pa- 
rameter U when the value of the Fermi level of the leads 
(i eq is taken equal to zero. We willingly choose to con- 
sider the particle- hole symmetric case (ed = —U/2) since 
we know that it is a delicate case in the sense that the 
EOM approaches developed so far have failed to describe 
it correctly. The density of states shows a three-peak 
structure as soon as U becomes larger than T, with two 
broad peaks and a narrow Kondo resonance peak. The 



C. Out of equilibrium 

1. Differential conductance 

Out of equilibrium, the density of states in the dot is 
greatly influenced by the bias voltage or the difference 
between the chemical potentials of the leads. Fig. [3] re- 
ports our results for the noncquilibrium density of states, 
again in the particle-hole symmetric case at T/T = 10~ 3 
and U/T = 4 for different values of the bias voltage V. In 
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FIG. 3: (Color online) Nonequilibrium density of states in the 
particle-hole symmetric case at U/T = 4 and T/T = 1CF 3 for 
different values of the bias voltage V . The chemical poten- 
tials of the two leads are taken equal to Hl/r = ±V/2. The 
Kondo resonance peak splits into two side peaks located at 
uj = ±V/2, i.e. at the positions of the left- and right-lead 
chemical potentials. For convenience we choose to represent 
the two energy scales (energy uj and bias voltage V) as nor- 
malized by the factor . 
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FIG. 4: (Color online) Differential conductance dl/dV ver- 
sus the bias voltage V in the particle-hole symmetric case 
at U/T = 4 for different values of temperature. The curves 
show a zero-bias peak, followed by the beginning of a broad 
Coulomb peak at large bias voltage. The differential conduc- 
tance is reduced when either the temperature or the bias volt- 
age increases, suggesting that the Kondo effect is suppressed 
by temperature or nonequilibrium effects. 



contrast with the equilibrium situation, the Kondo res- 
onance peak splits into two lower peaks pinned at the 
chemical potentials of the two leads. The reason is that 
the transitions between the ground state and the excited 
states of the dot are now mediated by the conduction 
electrons with energies lying close to the left- and right- 
lead chemical potentials. 



FIG. 5: (Color online) Differential conductance dl /dV versus 
the bias voltage V at T/Tk = 0.1 in the particle-hole symmet- 
ric case for different values of U. The inset shows that the dif- 
ferential conductance as a function of normalized bias voltage 
V/Tk scales to a single universal curve dl/dV = /(V/Tk)- 
At higher voltages, the universal behavior is destroyed by a 
broad peak resulting from charge fluctuations. 

We then compute the differential conductance as a 
function of bias voltage for different temperatures and 
plot the results in Fig. At low temperatures, the bias 
voltage dependence of the differential conductance shows 
a narrow peak at low bias (zero-bias anomaly) reflecting 
the Kondo effect (mind the logarithmic horizontal axis), 
followed by a Coulomb peak centered around the value 
of the dot level energy. Increasing temperature dimin- 
ishes the intensity of the zero-bias peak, meaning that 
the Kondo effect is destroyed by temperature. 

In order to discuss the universality of the dependence 
of the differential conductance on the bias voltage, we 
plot in Fig. [5] the results obtained at zero temperature 
for different values of the Coulomb interaction U. In the 
inset, the differential conductance is found to be a uni- 
versal function of the rcnormalizcd bias voltage V/Tk, 
independent of other energy scales such as U or T. This 
one-parameter scaling is obtained over a large range of 
V. Universality is lost around V > IOTk- Note that 
when V/Tk < 0.1, the unitary limit is not completely 
recovered for the differential conductance due to the nu- 
merical accuracy in the self-consistency treatment, as was 
already mentioned before. 

The physical origin of the destruction of the Kondo 
effect is the decoherence rates induced by the voltage- 
driven current. As we discussed in Scc lIII B"2l these ef- 
fects are well described by our EOM approach since it 
incorporates higher-order terms in t a . They originate 
physically from the energy-conserving processes in which 
one electron hops onto the dot from the higher chemical 
potential while another electron hops out to the lower 
chemical potential. Since the processes involve two elec- 
trons hopping in and out, the lowest-order contribution 
is fourth order in t a . These rates broaden and diminish 
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FIG. 6: Color plot of the differential conductance dl /dV as a 
function of bias voltage V and dot- level energy Ed for U/T = 4 
and T/T = 10 -3 . The contour of the Coulomb peaks delimits 
the Coulomb blockade diamond, separating areas with well- 
defined dot occupation number M ranging from 0, 1 to 2 at 
low V, and areas of charge fluctuations at high V. In the 
J\f = 1 central valley, dl /dV shows a zero- bias peak typical 
of the Kondo effect. 

the Kondo resonance peaks in the density of states as the 
bias voltage increases, see Fig. [3J Their effect leads to a 
decrease in the differential conductance when V > Tk, 
as was shown in Figures 0] and [S] We will analyse this in 
more detail later in this section. 

To further demonstrate that the method can work in 
a wide range of parameters, we report in FiglS] the dif- 
ferential conductance in the V — £d plane in a 3D-plot. 
The figure shows the usual Coulomb diamond defining 
inside the Coulomb blockade regime Af = 1, where N is 
the total occupation number in the dot (J\f = Y^<r n cr)- 
The boundaries of the Coulomb diamond are related to 
the values of the renormalized dot level energies ±Ed and 
±Ed + U (with some additional renormalization effects 
in the mixed valence regime). Within the Coulomb di- 
amond along the V — line, one can clearly see the 
zero-bias peak as discussed in Fig. [5] At zero temper- 
ature, the unitary limit 2e /h is almost reached at the 
particle-hole symmetric point (e<j = — 17/2). When the 
temperature is increased, the zero-bias differential con- 
ductance decreases at this point, leaving aside two broad 
Coulomb peaks corresponding to the alignment of the 
dot level energy with the chemical potentials in the leads 
(e d = -U and e d = 0). 

2. Comparison with other studies 

We compare our results for the differential conduc- 
tance with those obtained by other groups using time- 
dependent Numerical Renormalisation Group^i and an 
imaginary-time theory solved by using Quantum Monte 
Carlo 22 i 31 , and plot the results obtained for the bias volt- 
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FIG. 7: Comparison of the differential conductance dl/dV 
as a function of bias voltage V with the results obtained by 
AndersZi and Han£i for -2e d /r = U/T = 2.5 and T/T = 
0.008 < Tk/T (We are grateful to J.E. Han for providing 
us with his data points). Our curve is plotted for T = 
in order to compare the three results in the strong coupling 
regime. dl/dV at small V is slightly different in the EOM 
approach because its value for Tk is smaller. At high bias 
voltage, the results of the three approaches agree perfectly. 
A little unphysical bump is observed for the EOM result at 
V — U = 2.5F, when the chemical potentials of the leads 
are aligned with the resonant levels of the dot (px, = ea + U , 
Mfl = Ed)- 



age dependence of the differential conductance at zero 
temperature for comparison (Fig. [?])• One finds a quali- 
tative agreement at low bias voltages, when the system 
is in the strong coupling regime. In that regime, our 
method slightly underestimates dl/dV because it gives 
a smaller Kondo scale. The three curves join at higher 
bias voltages, where a quantitative agreement is found. 
A little local bump is observed for the EOM result at 
V = U, when the chemical potentials of the leads are 
aligned with the resonant levels of the dot (tt£ = e c i + U, 
= £d) ■ This is related to the fact that we used the bare 
D a (s) functions (|32|) in the non-Kondo regime, leading 
to divergence at ^l(r) = {£d,£<i + U}. This bump can 
be smeared out by introducing a finite width of order T 
into the D a (e) functions, as can be physically originated 
from charge fluctuations on the dot resonant levels. 



3. Crossover from strong coupling to weak coupling regime 

When a bias voltage is applied to the leads, it is inter- 
esting to know whether or not the decohcrence effects in- 
duced by the voltage-driven current (cf. Sec lIIIB"2j) may 
drive the system from strong to weak coupling regime. 
We point out that this problem has been discussed in pre- 
vious studies for the Kondo model using either a pertur- 
bative renormalization approach^ or a slave-boson tech- 
nique within non-crossing approximation^. We would 
like to tackle this question for the Anderson model with 
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two leads using the EOM scheme. 

At zero temperature and at F > T K , Q^{oo) for the 
Kondo regime behaves as 



In 



(w - V/2 + i^)(u + V/2 + i^) 



(,41) 



where 7^ is the decoherence rate induced by the bias 
voltage as given by Eq. piap . gi yen by Eq. (|4"Tj) 

develops a pole 1 ^ as soon as is smaller than a char- 
acteristic energy scale T* 



- V 2 /4 
1%/V 



V < V2T K 

V > V2T K 



(42) 



From there, we define a criterion controlling the crossover 
between strong coupling < T*) and weak cou- 

pling (7^) > T*) regime, as proposed in RefHol. In 
order to obtain the nontrivial decoherence rate 7^ 4 ) as 
a function of bias voltage V, we replace D a (e) of the 
bare J a a,0a' in the decoherence rate by the 'dressed' 
D a {e), which is identified with the denominator of the 
Green function (j!7[^ £. Thus we can define a renormalizcd 
Jaa^o-i = 2t aa tp a iD a {ii a ). In the limit V » T K , we find 

J a .,p*> oc l/[2 1h(V/Ta:)] and 7 M ex V/[2 ln(l//T K )] 2 , 
which is always larger than T* . The results for the renor- 
malized decoherence rate 7^/^if as a function of the 
bias voltage are reported in Fig. 8(a) for different val- 
ues of U . Strikingly, the curves for the different values 
of U coincide, underlining the universality of the evolu- 
tion of 7W /T K as a function of V/Tk- Combining the 
results for T* /T K and /T K , one can derive the uni- 
versal crossover bias voltage V c /Tk from strong to weak 
coupling regime. 

At finite temperatures, the derivation for T* is the 



same except for replacing 7 



(4) 



(4))2 



T 2 T 2 



Eq. (grj). The results are plotted in Fig. |8(b)| in the V—T 
plane, displaying the crossover from strong coupling to 
weak coupling regime. Although the physical mechanism 
at the origin of the crossover is different, both bias voltage 
and temperature drive the system to the weak coupling 
regime. 
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FIG. 8: (Color online) (a) Decoherence rate 7'*' and char- 
acteristic energy scale T* versus the normalized bias voltage 
V/Tk at T /Tk = 1CF 1 in the particle-hole symmetric case 
for several values of U/T. 7'*' /Tk is a universal function of 
V/Tk over a large range of V. The comparison of both en- 
ergy scales and T*) allows one to determine whether the 
system is in the strong coupling regime 

( 7 ( 4 ) < T*) or weak 

coupling regime (7" > T*). (b) Stability phase diagram of 
the strong coupling and weak coupling regimes in the V — T 
plane. The crossover tempeature T c (V)/Tk is a universal 
function of V/T K . 



4- Nonequilibrium occupation number in the dot 

Typically, at equilibrium and for zero temperature, 
(n a ) is mainly determined by the weight of the broad 
resonance peak far below the Fermi level. The narrow 
Kondo resonance near the Fermi energy has little weight 
in comparison. Thus, even if a EOM approach in a 
certain approximation scheme happens to describe only 
qualitatively Kondo physics, it is able to determine nu- 
merically the occupation number that agrees reasonably 
well with the Bethe ansatz or NRG. 

When the system is driven out of equilibrium, the 
problem becomes more complicated as one should use 



lesser Green functions instead of retarded ones to com- 
pute the expectation values. As discussed in Sec lIV Al 
the only place where this cannot be circumvented is pre- 
cisely for the dot occupation number (n a ) appearing in 
the Green function (|17[) . A rigorous treatment would re- 
quire to compute the lesser Green function Q< (oj) and 
then obtain (n a ) according to Eq. (|39| . which is beyond 
the scope of this work. As described in Sec lIV A| we 
used instead the Ng ansatz to compute the dot occu- 
pation number. On the other hand, if we consider the 
particle-hole symmetric case (e^ = —U/2), with a sym- 
metric bias voltage setting [^l,^r] = [V/2,— V/2], one 
obtains (n a ) = 1/2 by symmetry. However, we noticed 
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FIG. 9: (Color online) Occupation number in the dot (n a ) 
versus the bias voltage V in the particle-hole symmetric case 
for U/T and T/T = 8.51CF 4 , and under an asymmetric bias 
voltage setting — and /j,r = —V. 



that the calculation of the occupation number by apply- 
ing the Ng ansatz to our Green function leads to slight 
deviation from (n a ) = 1/2. In |Appcndix IV we show 
how to solve this problem. 

On the other hand, for an asymmetric bias voltage set- 
ting, the occupation number (n a ) is no longer fixed by 
symmetry arguments. Let us take \jj,l,Hr] = [0,— V], 
the bias voltage dependence of the occupation number is 
shown in Fig. [5J As V increases, (n a ) decreases rapidly 
till V passes the dot-level energy U /2 and comes to sta- 
bilize at large V. This can be qualitatively explained by 
the fact that at large V, the current through the dot no 
longer increases monotonously with the bias voltage and 
reaches a horizontal asymptote. This makes the occupa- 
tion number insensitive to the bias voltage. 



particle-hole symmetric case. When the dot is symmetri- 
cally coupled to the leads, the linear conductance reaches 
its unitary limit 2e 2 /h at zero temperature in the Kondo 
regime. 

We have also computed the nonequilibrium decoher- 
ence rate "/^ in the Kondo regime. At T ^> Tk, /Tk 
is found to be a universal increasing function of the nor- 
malized bias voltage V/Tk , depending on a single energy 
scale Tk- The scaling law holds over a wide range of V 
going from to IOOTr-. At low temperature, the density 
of states shows a splitting of the Kondo resonance into 
two peaks, pinned at the chemical potentials of the two 
leads. The height of the two peaks is controlled by the 
decoherence rate. 

As far as the differential conductance is concerned, it 
shows a zero-bias peak at low temperature, followed by 
a broad Coulomb peak at larger bias voltage. At low 
bias voltage, the differential conductance also obeys a 
universal scaling law as a function of V/Tk- Finally we 
have discussed the role played by the decoherence rate 
7 < - 4 ^ in driving the system from the strong coupling to 
the weak coupling regime. We have derived the crossover 
line T C (V) separating the strong coupling regime to the 
weak coupling regime. 
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V. CONCLUSIONS 



We have presented a study of the nonequilibrium ef- 
fects in the two-lead Anderson model. The calculations 
have been performed within a self-consistent EOM ap- 
proach generalized to the nonequilibrium situation. The 
approximation scheme presented in this paper goes be- 
yond the previous truncations of the equations of motion 
done at the second or fourth order in tunneling t a , by in- 
cluding contributions from the next orders (sixth order), 
which have been shown to be of great importance out of 
equilibrium. 

The situation at equilibrium is used as a benchmark for 
the approximation. The results for the density of states 
and the linear conductance at equilibrium are found to 
be quantitatively improved compared to those obtained 
by the EOM method using the Lacroix approximation. 
In the Kondo regime for instance, the Kondo temper- 
ature Tk is closer to the exact results found with the 
Bethe ansatz and NRG, and non longer vanishes in the 
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Appendix I. DERIVATION OF EQUATIONS OF 
MOTION FOR FINITE COULOMB 
INTERACTION 

A. 

In Scctionlnl we have derived the first equations of mo- 
tion. In this appendix, we present the detailed derivation 



of the higher hierarchy of equations and the decoupling 
scheme that follows. We derive the EOM of the higher 
Green functions on the right-hand side of Eq. J5J by using 
Eq. (g|). They are 



u-.k((n 9 c k(7 }) 

U S ;ak{{ftcksU)) 
(u>k:a 9 ~ U) {(c\-f 9 f a )) 



ta((n<rf*)) + ^ts\{{flck'sCka)) - {{c\, s f <jC k<J )) 
k' 

UUka) +t 9 ((n 9 f a )) +^2[ta((ftc k ^C k/(T )) - t s {{c\, s C kS fa)) 
k' 

{c\-J 9 ) - + E [**<<4 W«r>) + M<4/^>) 



(A.la) 
(A.lb) 
(A.lc) 



where we denote for a shorthand 

U a p...:ab- = w + e a + ep H e a — Sb — ■ ■ ■ , 

with {a/3 ■ ■ ■ , ab ■ ■ ■ } being any set of parameters within 
k's and cr's. 

For most practical purposes, the truncation is per- 
formed at this level by decoupling the second-order terms 
on the right-hand side of Eqs. (|A.1[) . see the paragraph in 
the main text after Eqs. (fTU|) . This is done by grouping 
all possible same-spin pairs of lead (c) and dot (/) elec- 
tron operators since we assume the spin quantum number 
is preserved through tunneling: any correlation between 



electrons of different spins has to come via the Coulomb 
interaction. A solution obtained at this level by neglect- 
ing the connected Green functions is exact to second or- 
der in hybridization^ 8 .. Numerous such solutions can be 
found in the literatu re 25 ' 26 ' 38 ' 39 , with some more elabo- 
rate than the others. 

However, as discussed in the main text, stopping the 
flow at this point will raise terms suffering from loga- 
rithmic divergences. In the following, we show how to 
go beyond the second-order to derive higher equations of 
motion exact up to the fourth order. To begin with, we 
consider the following second-generation EOM: 



Ug; kk >{{ftck>gCko)) 
^k':Sk{{c\, s f S C k a)) 



-U{{n a fic k , & c ka )) +ta{{n 9Cka )) +t a ((fic V9 U)) - ^ M(4"fr c *'» c *<r»« (A.2a) 

fc" 

ft k + U{{n s <b s Ck*U))+t s {{<b s f 9 U))-t^ 

fc" 

U((n a cl s f 9 c ka )) - t 9 ((n 9 c ka J^)) + t a ((cl,J 9 f a )) + Y^M(ct 9 Ck'' 9 c ka )). (A.2c) 



where we denote f%, k = {c} k , a Ck<j)- 

We proceed to insert the above Eqs. (|A.2[) into 
Eqs. (|A.lj) . The right-hand side of the latter equations 
involve new Green functions generated via the Coulomb 
interaction and others of the same hierarchy as the left- 
hand side. The latter Green functions can either move to 
the left-hand side or vanish in the wide-band limit since 



upon summing over fc, all denominators have poles in 
the upper half complex plane. Furthermore, we decou- 
ple {{c\c s c a )) w (c 9 c 9 )({c a )) and then use Eq. ((6]), since 
this decoupling should be exact up to order of t\ and for 
another reason which will be clear later. We end up with 
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[cj. ft - E la (uj :k )] {{n s c ka )) 



ta((naD) - U^tg ]^^ kk , {{n a fl c k ,„c ka )) + u k ,] 9k {{n a cl,-fec k(T )i A,.3a) 

k' 

fe' 



-17 



fe' 

+ C/ E [*<TWfc : fffc'(("^ c L/ff c fc'^» + t *Uk- 1 ak>(( n * C k9 C k'<rf<T)) ■ ( A -3c) 



where It is interesting to notice that at this level the prefactor 

of the Green functions on the left-hand side acquires non- 

£» = (A.4) interacting self-energy terms, while new Green functions 

k remain on the right-hand side. This allows us to focus on 

Sicr(w k ) = ^""^ t- [oJ~ kk / +^ k i 1 - k ] (A 5) the new Green functions, which is of most importance. 

4^ <T ' We list below these third-generation equations of motion, 

^Mm) = E ft"*** + - (A-6) 
fe' 

S 3CT (^:) = E [4^ife' + 4^:k'] ' (A.7) 



(Wff:fc'fc + U) {{n a flc ks C k 'a)) = -{.flc ks flc k , a ) + tg {{n a Tl s C kl a )) - ({c\„ g f a fffC ks C k > a)) 

fe" 

+ ^2 t <?{{.fl c k"<7.flcksCk>v)) -}Jg((n a cl„^c k9 c k 'g}}, 
fe" fe" 

(wfc' :ff fe - J7) {{n a c\,-f s c ka )) = -(c{, s f s flc ka ) - t s ((n a n s c ka )) -^t a ((cl„ a f a c\, s f s c k<J )) 

fe" 

+ ^to((fic k „ a cl, s f s c k<J )) + ^t s ((n a cl, g c k „„c ka )), 
fe" fe" 

(w fe ' :CTfc - U) ((n s cl, s c ks f a )) = (n s c\, s c ks ) - ^ M(4"*/fr4'er<W<>->) 

fe" 



(A.8a) 



(A.8b) 



(A.8c) 
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Up to order t*, we decouple Eqs. ()A.8|) by considering the following decouplings: 



ifa c kafa c k'a) 

( c fc'g-/ff /<r c fco") 
((ci"<jf<?fl C k5Ck><j)) 
{{flck"i7ftckgCk><j)) 
((riaC^-CkaCk'a)) 

( ( C k" a f" °k' S fa c ka ) ) 
{{flck"acl, s f S C ka )) 

{{ C l"sfsc\'S C kS.fa)) 

{{fUl'S C k"sC k sfa)) 



{fic k - a ){flc k ,„), (A.9a) 

{n*)m k -{fU k& )(cl,-U), (A.9b) 

(c{ rff h){ftc ka ), (A.9c) 

-fZ»k>{{fUk*U)) - (fU^dcl^U)) + f%„ k ,(fUk*)((U)), (A.9d) 

(/|c fcff )((/t Cfc , v c fc v)), (A.9e) 

/f„ fc ((n CT c fc v», (A.9f) 

-fv>k((<i*Uf*)) - (ol & U){(cl a c ka U)) + m, k {c{^U)((f a }}, (A.9g) 

{cl, s UWUk»„c ka )), (A.9h) 

-fZ» k ((claUU)) + flM'sUU)), (A.9i) 

flk»{{fUk & f*)) ~ flkfJjlcv.sU)), (A.9j) 

/£ fc «n»<v ff » - (/ic te )«4 ff / ff cvv)). (A.9k) 



Note that the last terms in Eqs. (|A.9dl |A.9g[ ) are added in order to avoid double counting from the first two terms. 
On the other hand, since {{n s n a c ka )) vanish in the wide-band limit after summing over k, they are removed from 
now on without further notice. Eqs. (|A.8|) then become 



{^s:k'k + U){{n a flc kg c kl<7 )) = ^Uf^vdflcksU)) + u :k ,(flc k s) -^tsfw-k {{n a c k 'a)) 

k" k" 

-{fUks^^fk-k'df"))' ( A - 10a ) 
k" 

{Uk>:ak-U)(KcUfaCka)) = E ((4>M*)) + (4'sU) + E WfcV] 

k" k" 

-<4'*/*> (A.lOb) 
k" 

(w*:«r* - U){{n s cl, s C ks U)) = -{fUks) [{c\,sU) + E + E ^'M S UU)) 

k" k" 

-£WSk»<</£W«r>) + {Uu, - U).f k , k {{n s U)). (A.lOc) 



To obtain Eqs. (jA.lOal IA.10b[) in a more compact form, we have used the following equation of motion, 



>:k((n a c ka )) = - Jlc ka ) +^2t a [((cl, a c ka f a )) + {{fic k : a c ka )) 



(A.11) 



18 



and taken Eq. ([8]) for Eq. (|A.10c|) . Introducing Eqs. (IA.10[) into Eqs. (|A.3|) and after some straightforward algebra, 
we eventually obtain 



[u !k -S lo {u :h )]((n 9 c kir )) = t a ((n s f a )) - E t a t s [D S:kk , + D k , 9k (cl 9 f s )\ fZ„ k ((f a )) 

k'k" 

+Y, 5rf (uj. k )((n rT c ka )) + E t a t s D s .. kk , fk"k((fs c k'sf a )) 

k'k" 

+ t<,UD k ,. &k fZ„ k ((c{,^f a )), (A.12a) 

k'k" 

[iO S:ka - £ 2<r (u, :fe )] ({flcteU)) = + (fUka) [l -^UD^k^gU) +X)*"«4'a/ffCfc"«r)>) 

fe' fc" 

- E + E £AW£»tf] ((A)) + E 4Dk':,kf^ k ({cUfafa)) 

-J2 t ° D ^k'[j2 ts -fk"k - U:k>{fUke)]((n<7C k , a }), (A.12b) 
fc' fe" 

[ Wfe:aff - 17 - £ 3CT (u; fc: )] ((cUUf*)) = -t s ((n s f a )) + {c\ s f s ) [l + J2^D k:akl (</Jc fc ^) + E «»{(/W»)>) 

fe' fe" 

+E + (cL/^E^^^'/fe-fe'] + E tiD k .. ak 'f kk „{{fU k ' S f a )) 



[w-k - Si5-(o; : fe)] ((n a c krT )) 



-E^fc:fffc'[E^^" + W ^'( C L/ff) ((™<x C fc'<r)), 

fc' fc" 

-(/>fc CT > + E kA*** [EWa - w,* (/]c feCT >] ((n ffCfc / ff )) 
fc' fc" 

+S 5CT M((n^Cfc CT ))+E^/fc'fc((/<r»- 



where 



D a p—:ab— — —Uuj a p^. ab ^(uj a p... :a f } ...±U) , 
The sign in front of £/ is the same as the sign in front of e a in uj a p... :a b.... The self-energy corrections are 

E lff (w :fc ) = E^KL' ' 



fe':<rfej ' 



E 2(T (w :fc ) EE E 2ff ( W:fc ) + E [4^:fcfc'/fc"fc'-4^fc'-fc/fc ff 'fc" 
fc'fc" 

E3a(Wfc) = S 3 «tH:)-E [^fc:5fc'/fc"fc'+4^fc-fc'/fc"fc' 



(A.12c) 

(A.12d) 

(A.13) 

(A.14) 
(A.15) 
(A.16) 



We also define two following functions, 

£ 5 <r(W:fc) = E t i[ D k'-<rkfk'k" ~ Ar:fcfc' /fc"A 



k'k'' 



+ u, k Y,t*[D k ':„k{c{,J a )+D a:kk '{flc k ' rJ )\ 1 . (A.17) 



After truncation there appears a new Green function ((n a c ktT )) on the right-hand side of Eqs. (|A.12allA.12cj) . which 
has to be calculated separately. Its equation of motion, given by Eq. (|A.12d[) , is derived in | Appendix I B| Note that 
in deriving Eqs. (jA.12bMA.12d)) . we approximate 



(A.18) 



In doing so, we assume that the lead electron energies s k , e k i ~ Hl(r) cancel each other; this assumption is seconded 
by the numerator f^ k , = 5 kk i fp{s k ) at zeroth order. Thus it is valid to use Eq. ([7]). This approximation should not 
affect the density of states around the Fermi level. 

We emphasize that Eqs. (|A.12[) shown above are exact up to fourth order. The prefactor of these equations of 
motion acquires second-order corrections. Their imaginary part takes different values in the different regimes of the 
Anderson model, as is discussed in Scc lIIII 
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B. Derivation of the equation of motion of {{n a Cka}) 



Here we expand the equation of motion of {(n a c ka )). The equation is already given by Eq. (| A. 1 1|) . We now derive 
the next higher equations 

^a:kk'((f^Ck'aCka)) = ~ ^ t a ( (cj^Cfc^ C ka ) ) + t a ( (n a C ka ) ) - t a ( (n a C k ' a ) ) - U( {n s /* C k ' a C ka ) ) , (A. 19a) 

k" 

Uk':ak((c{, a c ka f a )) = f k , k + t a {{n a c ka )) +^t a {{c\, a c ka c k ^ a )) + U{{n s c\, a c ka f a )). (A. 19b) 

fe" 

Similarly at this stage, we go on to decouple ((4" CT c fevCfe<r)) fZny{{cka)) — f k " k (( c k'a))- After removing terms that 
will vanish in the wide-band limit upon summing over k, Eq. (|A.11|) becomes 

[w iJk -S lff ( W:fc )]((n ff c fcff )) = -(/K^ + E^-^^'fc + 

it' 

+ [wfe,U((^4' < r C ^/<r>)-^L'(K/i C fcVC fcCT >)] • (A.20) 

fe' 

where SiaC^fe) = /J^^CT-fefe' w fc~'-o*] ■ Next we derive the equations of motion of the Green functions on the 

fe' 

right-hand side of Eq. (|A.20j) 

(oJa-.kk' +U){{n s flc k , a C ka )) = -^t 9 ((& k „ s Uflck' a Ck a )) + ^2ta((Uc k "afaCk'aCka)) 

k" fe" 

+t rT ((n g -n a c ka )) - t,j{{n s n a c k > a )) - t a {{nac\„ ^c^ a c ka )) , (A. 21a) 

k" 

{^k'-.ak - U){{n s c\, a C k<J f a )) = (n s cl, a C ka ) ~ y^^((4"g/^ C fc'<7 C fc^/^)) + *cr( vWCfctr)) 

fe" 

+ X^</cUfe»a4'<7 C W<r>> +X>^ C fc'a C ^ C >^}}- (A.21b) 



We now decouple Eqs. (|A.21j) according to the following decouplings: 

(n s c\, a c krT ) w fk'ki n s)> 

((cl"&f5fl C k'<yCka)) ~ (/Jc fc / (T )((4„ & / & C fc(T )) - (/Jc fce r>((4„ e ./ S Cfc/ er )), 

((ftck"aftck'aCka)) ~ (/aCfev)((/ffCfe"aCfe CT )) - (f^C k a}((ftc k "aCk'a)), 

((fifiv*^)) ~ fk"k'(( n » c ka)) - fk"k(( n sCk'a)), 

{{ c l"sfsA'a C k<yf<y)) ~ .fk'k{{ c l" S fsfa)) ~ (C k ,ja}{{cl„ & fcc k( ,}}, 

((/ffCfe"a4'a C fc<^/^)) ~ fk'MfUk"sfa)) ~ {c\, a fa) {{fUk" sC ka )) . 



(A.22a) 
(A.22b) 
(A.22c) 
(A.22d) 
(A.22e) 
(A.22f) 



Again the Green functions ({ng-n a c k(T }) vanish in the wide-band limit and are removed hereafter. Using Eqs. (|51 IA.lal) . 
Eqs. (|A.21|) become 

(w^fe' +U){{nvflc k > a c ktT )) = [^ia/fe"fe-w :fe '(/]c feCT )]((naCfev}}- [XX/k"fc' ~ u-.kUUk',,) ((n ff c fe( ($,23a) 

fe" fe" 

(Wfe': CT fe - C/)((^4 vCfeCT / ff )> = tf, fc ( W:<y - U){(n S U)) - [Y,t°fk>k» +>":k{<bM] ««*<*„». (A.23b) 



Combining Eqs. ([A~20l 1X23)) and using Eq. (jA~T8)) yield Eq. (|A.12d|) . 



Appendix II. DERIVATION OF EXPECTATION VALUES 

At equilibrium, the hermiticity of expectation values holds, e.g. (Cf. /(T c ka ) = (4o- c ' ! ' CT )' (fa c ka) = { c \ a f<?) because 
= — /f(w)[£/£(w) — ^(w)]. In the Keldysh formalism^. Gf(w) is the lesser Green function, while Q^ a \ui) is 
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the retarded (advanced) Green function. This relation is nothing but the spectral theorem or equivalently Eq. (|38l) . 
with which it is standard to transform expectation values into a functional of the retarded or advanced dot Green 
function^. However, the spectral theorem does not apply out of equilibrium^ and it is therefore necessary to invoke 
the nonequilibrium Keldysh formalism^. We show how to rewrite the expectation values in Eq. (|17|) in terms of 
integral functions. We calculate for our purposes the following expectation values: 



f duj „^ , , 



-itn 



f F {e k )g a a {e k ) ■ 

M = 



2^ 



2ni tu-k + iS 



dio 



f F {e k )g r Ae k ) 
du> „^ 



it 

duj g<{u) 



2ni ui-k — iS 



(B.l) 



(B.2) 



-fe'er 



Cfco 



2jt 

doj f , 



( W ) 



+ 3 LHS<Hg fe Q <» + fl£>)0>)<£«>)] } 
r/F(e fe )^(e fe )-/f(e fe ')^(^') 



dio 



5<H 



(B.3) 



where flll (o>) 



± £5) 1 , g ka {w) = 27rz/p'(efe)(5(aJ:fe) are the bare lead Green functions. Using Eqs. ([B.HIB.3I) . 
we obtain after summation two related functions defined by P a {oj) and Qo-(lu) 



P,(w) 



E 

fe 

E 

fefc' 



Ckn 



£ k + iS 



E ■ 

a=L,R 

E : 

a=L,R 



r, 



iS 



otu 
7T 



ir„0?(e)] 



(B.4) 



(B.5) 



Notice that the imaginary part of the denominator is always positive, so that the pole e = uj + iS remains in the 
upper half complex plane. In deriving Eqs. (|B.4IIB.5|) . some terms, particularly those associated with g^(cj), vanish 
in the wide-band limit, since upon summing over k, all denominators have poles in the upper half complex plane. 
Hence we have shown that in the wide-band limit, the nonequilibrium functions P a {ui), Q a (uj) take the same forms as 
in equilibrium, except that the left and right leads have different chemical potentials. No knowledge of lesser Green 
functions is needed. This constitutes a huge simplification in the computations. 



Appendix III. CHARGE CONJUGATION SYMMETRY 

This appendix is devoted to proving charge conjugation symmetry of the dot Green function given by Eq. (|17l) . 
We follow the scheme established by V. Kashcheyevs et who proved that this identity holds for the Lacroix 
approximation. The Anderson Hamiltonian ([1]) attains its original structure if replacing the particle operators by the 
hole ones, /j = f„, c\ a = c ka , along with 

C{e a ) = -e a - U, C(U) = U, C{t a ) = -t*, C(s k ) = -e k , C({n a )) = 1 - (n CT ), (C.l) 

where C is the charge conjugation operator transforming electrons quantities into hole ones and reversely. The hole 
dot Green function is related to the particle dot Green function by charge conjugation symmetry 

c km] = ((aj a ),aji))) u = -g^-ui). (c.2) 

Eq. (fl7|) obeys this symmetry if the following rules are respected: 

C [ui a (uj)} = U - Uio-(— w), C [u 2a {uj)] = 1 - U2o-(-w), 

C [£° («)] = -E° (-w), C[P <T (u)] = -P v (- U ), (G3) 

C[Q a (w)] = Q„(-u) 
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Therefore we deduce that C [P a (u> a:r7 )} = -P CT (-aw), C [Q a {u> w)] = Qo-(-kw) - £° (-ow), and similarly 
C [P a (-LJ :cr cr + U)] = -P a (oj aa: +U), C [Q a (-oj. aS + U)] = Q a (oj a9 . + U) — £° + U) . For the Lacroix approxima- 
tion, this is enough to prove that the Green function respects charge conjugation symmetry C [£/£(w)] = — Q a {— w)^. 
However, Eq. (jTTJ) has additional self-energy corrections in the arguments of -Po-(w) and Q a (uj). Here we will show 
Eq. jG3} still holds for Eq. ^ . Using Eq. del}, we find 



C [Ua:kk'] = U k k':<x ~ U, 



C [uJk':ak] = ^ak:k' + U. 
C [Dk:ak'\ = D k:ak 1 



UJ — > — UJ 7 



c 



S 2 (w : fe) 



-E2o-(— Wfc:), 



S 3 (w fe: ) 



Thus if we redefine 



L sJk'k 



kk , ^>S:ka ~ ^2a{^>:k) 

p-oW.o) = 2^ ~— — -, 

k ^S:ka ~ ^2a{^J-.k) 



Qa{-W.. aS +U) = J2 
kk' 



-T, 3lT (-LJ :k ). 



i sJkk' 



kk , w k .. aS + U + £ 3cr (w fc: ) ' 
-Wk-.od + U + £ 3cr (w fc: ) 



(C.4) 



(C.5) 



One can see that the above functions obey the above transformation rules ()C.3j) by analysing the transformation rules 
of their self-energies T,^. In Eq. ([TT]). there are yet another two terms, defined by 



kk 



■f [w : fc-£6 CT (w : fc)][w.fc-£i s .(w : fc)] : 



t a ^5a(0J. k )(f^C kcr ) 



[uj :k - Y, 6a (uj. k )}[uj. k - Ei ff (w :fe )]' 



(C.6) 



kk' 



One can show that 

C[P la {w)]=-P la {-w), 

C [Qla(u)] = Qla(-U) ~ £° (-W) 

in the wide-band limit. With this assumption, performing Eq. (|C.ip on Eqs. (|12ll5p yields 

C [E 1<T (U; k )] = -Ei ff (-w fc: ), C [S 5o -(cj :fe )] = -£ 5cr (-cj fc: ), C \E 6 (uj : k)] = -£ 6 (-w fc: ). 

Therefore Eq. (|C.3I) holds for (5i CT (w) and Pi a (u>); similarly for u\ a {ui) and U2a{w). As a result, we prove that Eq. (fTT|) 

maintains charge conjugation symmetry by obeying the particle- hole relation Eq. (IC.2I) . 

| 



Appendix IV. SYMMETRY WITH RESPECT TO 
THE PARTICLE-HOLE SYMMETRIC POINT 



In the previous section, we showed that the Green func- 
tions calculated from electron and hole Hamiltonians are 



related by charge conjugation symmetry. By changing 
f a — > C{f a ), Cka — > ^cjL), the Anderson Hamiltonian 
for holes C [H] = Hh is, within a constant energy, 



J 



Hh = J2(- £ ^ C ( c lka)C(.c ak<T ) + J2 (-e«r - U) C(n a ) + UC(n t )C( ni ) - £ (t aa C(cl ka )C{f a ) + H.c) , (E.l) 



where we keep the parameters of the original Hamiltonian (|E.2[) of another system, whose parameters share with 
for electrons. It maintains the structure of an Anderson Hi- 
Hamiltonian, but with transformed hole operators. 

Here lies another symmetry: to the hole Hamiltonian 
Hh (|E.1[) corresponds an electron Hamiltonian Hf yst2 



n s j st 



U) n a + Un tni - J2 (t^ikJa + H - C 



(E.2) 
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System 1 



U 



R 



System 2 



L e d 



FIG. 10: Schematic representation of two dual systems in 
the electron picture. The dotted line represents the particle- 
hole symmetric point E* . Each system is equal to the mirror 
symmetry of its dual counterpart around E* in energy space 
. The electrons in System 1 must behave exactly the same as 
the holes in System 2, and reversely. 



We call dual systems two systems showing the symmetry 
^Systi = ySysti^ as snowri i n Yig. rjUJ For instance, one 
can have the following parameters 





System 1 


System 2 


He 


£d 


= -2 ; U = 6 


e d = -4 ; U = 6 


n h 


C(e d ) 


= -4 ; C(U) = 6 


C{e d ) = -2 ; C(U) = 6 



An electron in the first system behaves exactly as a hole 
in the second (dual) system and reversely. 

This symmetry is slightly broken by our approxima- 
tion scheme; the worst case is in the Af = 1 CB regime. 
The reason could be due to the fact that at order tt we 



do not treat the particle and hole contributions on an 
equal footing. Therefore, because the transition rates of 
the self-energies £20- and £3^ have different values in the 
frequency range e a < uj < e a + U, it leads to slightly 
asymmetric renormalization and broadening of the res- 
onant peaks at e a and e a + U, as can be shown in the 
particle-hole symmetric case, which affects the occupa- 
tion number. For instance, at the particle-hole symmet- 
ric point (e a = —U/2, = — /Ar) ; the dot occupation 
number (n a ) is expected to be exactly 1/2 in equilibrium 
or in the symmetric bias setting. Our numerical result 
shows deviation by a few percents at worst. However, 
it has almost no effect on the low-frequency density of 
states structure. 

In order to restore the symmetry, one compute the 
Green function in the dual system. Because of the defi- 
nition of the duality, we have the identity 



(E.3) 



where Systems 1 and 2 are dual of each other. Using 
charge conjugation symmetry (cf |Appcndix III ) on Sys- 
tem 2, we can express this equality in terms of electron 
Green functions only, that is 



->Systl 



(w) = - [g* yst2 (-u;)Y 



As mentioned earlier, this equality is slightly violated 
at high frequencies by our approximation scheme. We 
therefore symmetrise the two by setting 

0» = - [^ yst2 (-^)f }/2. (E.4) 
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